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Abstract. We establish sharp L^-Sobolev estimates for classes of pseudodiff- 
erential operators with singular symbols 1291 1221 whose non-pseudodifferential 
(Fourier integral operator) parts exhibit two-sided fold singularities. The op- 
erators considered include both singular integral operators along curves in 
with simple inflection points and normal operators arising in linearized seismic 
imaging in the presence of fold caustics 1321 1101 [TT] . 



1. Introduction 

This paper is concerned with L^-Sobolev estimates for some classes of operators 
sharing a common microlocal geometry that displays two different kinds of sin- 
gularities. One of our results is that in certain cases these singularities interact, 
producing worse estimates than one might otherwise expect. We begin by de- 
scribing the main features of this geometry and the motivations for studying these 
operators. 

1.1. The cubic model. In ]R^, consider the Hilbert transform along a model curve 
j{t) :— {t,t^) with a simple inflection point, 

(1.1) Hf{x)=p.v. f f{xi-t,X2-t^)^. 

It is well known that H (and similar but much more general operators) are bounded 
on L^(M^) (see, e.g., |45l |46l E].) The Schwartz kernel of "H is given by 

(1.2) Ku{x,y) = 5{x2 - y2 - {xi - yif) x p.v. ( — ^ — ) , 

which incorporates both pseudodifferential and Radon transform type singularities. 
In fact, the wave front set of K-^ satisfies the inclusion WF{K-h) C A' U Cq, where 

• A = {(x, i]x,£,):x& M?, {0}} is the diagonal of T*M? x T*M? , 

• Cq is the conormal bundle of {x2 — 2/2 — {xi — yi)^ — 0}, and 

• ' denotes the usual twist map in microlocal analysis, sending 
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The same geometry is present for fractional integral operators along 7, 

f{xi-t,X2-t3)j^, 0<a<l. 

Similar operators, combining fractional integral singularities along the curve with 
fractional derivatives transverse to it, have been studied for various curves in 
or in nilpotent groups [HI [TU |3] . Operators with kernels that display two different 
types of singularities arise in a variety of settings and go under various names: In 
the terminology of |36j . H is a singular Radon transform, while in the language 
of |29l I22j , and JTq are pseudodijferential operators with singular symbols. The 
latter concept was originally introduced to describe parametrices for variable coef- 
ficient wave and other real principal type operators, but has since found numerous 
applications in integral geometry [THl HH] and inverse scattering [201 US] • 

In the present paper we obtain sharp L^-Sobolev estimates for a class of operators 
containing H and J7q, and with similar microfocal geometry. One feature is that 
in some cases the estimates are worse than one would expect by considering the 
orders of the pseudodifferential and Fourier integral operator parts of the operators 
separately. For example, Ja '■ -ffcomp — ^ ^100"^ ~ ^y^' rather than for 

r = max(a — 1,— |). Such behavior can already be seen, for general families of 
curves (for a close to 1) in Greenblatt [T5] . 

However, the Schwartz kernels of the operators we are most interested in are not 
associated with a nested pair of submanifolds [18] in the same way that those of 
fractional integrals along submanifolds are. We thus need to formulate this class of 
operators using a microfocal approach; from this point of view, 
and JTq, e '"~2 (A, Co), where, for p,l e M, the space /P''(A,Co) denotes the 
class of operators on whose Schwartz kernels belong to the class of paired 

Lagrangian distributions associated with the cleanly intersecting Lagrangian man- 
ifolds A' and Co' in T*{R^ x R% We refer the reader to [H HI [Ml HH] for the 
theory of the 1^'' classes and to |2]for the basic definitions that we need here. The 
key facts relevant for this paper are that 

(i) The canonical relations A, Cq intersect cleanly in codimension one [18] . 

(ii) Co is a folding canonical relation in the sense of Melrose and Taylor [5S] or a 



two-sided fold in the sense of [Tni[T7] (see §^2.1 for details) 



(iii) The fold surface of Co equals its intersection with the diagonal, A n Co. 

We consider the pair (A, Co) mainly as a stepping-stone to the seismic geometries 
described below, and do not study non-translation invariant versions of H and J^a 
along the lines of [HI [T4l|42l|4]. 



1.2. Single-source Seismic Imaging. A very similar geometry occurs in a com- 
pletely unrelated problem, involving forward scattering maps in linearized seismic 
imaging [351 [TUl [H] . Let Y = := {1/3 > 0} model the subsurface of the earth, 
and dY = {j/3 = 0} its surface or the surface of the ocean. Let E C dY x dY denote 
a source-receiver manifold of pairs (s,r), and T > the total time-length of the 
seismic experiment. The data space for a typical such imaging experiment is given 
by the collection X — T, x {0,T), which corresponds to surface measurements being 
made at locations r and times t, < t < T, resulting from seismic events with 
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(idealized) delta-function impulses at sources s and time to — 0. The dimension 
of the data space is of course one more than that of the source-receiver manifold. 
It is known that under fairly general conditions, the corresponding (formally) lin- 
earized scattering operator F, mapping perturbations in the sound speed c{y) in 
Y (about a known smooth background co(y)) to perturbations in the pressure field 
measured on X, is a Fourier integral operator (hereafter referred to as FIO) as- 
sociated with a canonical relation C C (T*X \ 0) x {T*Y\0) gOl US IM] • More 
specifically, for some order m S K dictated by the dimensionality of the problem 
at hand, F belongs to the class of FIOs, I'^{X, Y; C), consisting of operators map- 
ping £'{Y) — > I?'(X) whose Schwartz kernels are Fourier integral distributions 
associated with C" cT*{X xY)\ 0. 

The L^-Sobolev regularity of the operator F, which in turn is tied to the question 
of microlocal invertibility, turns out to be heavily dependent on the data collection 
geometry and certain non-degeneracy hypotheses regarding the background sound 
speed coiy), assumed smooth and known. Historically, the first situation to be 
mathematically analyzed was the single source geometry without caustics, for which 
(a) the seismic data is obtained from a single source, Sq £ dY and the receivers 
range over an open subset of dY (i.e., dim X — dim Y — ?>)\ and (b) the rays 
flowing out from the source sq have no caustics. The latter assumption means the 
following. Let us consider the smooth conic Lagrangian A^^ C T*^ \ obtained 
from the flowout of the points (^Oj^o) in the characteristic variety of the wave 
equation above the point (x,i) = (sqjO) using the Hamiltonian vector field of the 
principal symbol of the background wave operator, Dco- There will be no caustics 
if the spatial projection from this Lagrangian manifold. Try : As^ — > Y , has 
full rank everywhere (except of course (so,0)). Beylkin [3 showed that under 
these assumptions C is a local canonical graph, so that F*F is an elliptic pseudo- 
differential operator (hereafter referred to as vJ/DO) on Y . Application of a left 
parametrix of F*F then results in high-frequency linearized seismic inversion. This 
conclusion also holds for other data geometries, under the traveltime injectivity 
condition [251 IMl which ensures that the canonical relation C satisfies the 
Bolker condition f3T|, i.e., ttl '■ C — > T*X \ is an injective immersion; this was 
further weakened in some situations pTj- 

However, caustics are unavoidable in physically realistic velocity models 33" and, 
since they typically lead to artifacts in images, understanding the structure of F and 
F*F when the background sound speed exhibits caustics is a fundamental problem 
in exploration seismology. Nolan [32] showed that, for the single source geometry 
with caustics of at most fold type, C is a folding canonical relation. It was then 
shown that, for any folding relation Cfoid and any F £ I"^(Cfoid)j 

(1.4) /2"'0(A,Ci), 

with Ci having the properties (i,ii,iii) as in §^jri] above ([32], Felea [TU]). In partic- 
ular, for the seismic problem, F*F e P'^{A,Ci) . 

1.3. Marine seismic imaging. Another seismic data set of interest comes from 
marine (or offshore) imaging. A mathematical idealization of the experimental 
setup is as follows. A survey vessel trails behind it a cable containing both an 
acoustic source and a line of recording instruments. The point source consists of 
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an airgun which sends acoustic waves through the ocean to the subsurface. Re- 
flections occur when the sound waves encounter singularities in the soundspeed in 
the subsurface, such as discontinuities at interfaces of sedimentary layers, and the 
reflected rays are then received by a linear array of hydrophones towed behind the 
vessel. The vessel then makes repeated passes along parallel paths (say, parallel 
to the xi axis) contained in an open U C dY. In other words, the source-receiver 
pairs form an open subset of 

S = {(r, s) <eU xU -.r ^ (ri,r2,0),s = (si,S2,0),r2 = S2}, 

which is a codimension one in dY x dY. Thus the data set is overdetermined, with 
dim X ^ A> dim F = 3. The forward operator F e I^^^{X, Y; C), and in [TT] 
the first two authors identified the structure of C, under a natural extension of the 
fold caustic assumption to this context. General canonical relations having this 
structure, along with some additional nondegeneracy conditions, were called folded 
cross caps, and a composition result was proven: if C is a folded cross cap and 
/™-3(X,y;C), then 

(1.5) ^*Ae /2™-5.5(A,C2), 



with C2 C (T*Y \ 0) X {T*Y \ 0) again satisfying (i,ii,iii) from §[1.1 In particular, 
for the marine seismic imaging problem, F*F e /2'2(A,C2). 

Now we can compare the normal operators for the single source and the marine 
geometries (both in the presence of fold caustics). It follows from [29j [22| that 
microlocally away from A n C2 and for any p, Z € M, 

(1.6) P''(A,C2) -^F+'(A\C2) + F(C2\A). 



From (1.5) and (1.6), one observes that since I — ^ ior the marine geometry, the 



order of the non-^DO part, i.e., the FIO part, of the normal operator F*F is 



lower than the ^fDO part. On the other hand, from (1.4) and the analogue of (1.6) 
one sees that for the single source geometry the two orders are the same. In other 
words, the artifact resulting from the FIO part of the normal operator for the single 
source data set is at least as strong as the ^DO part and hence is nonremovable, 
but the fact that the artifact for the marine data set is | lower order (away from 
A n C2) leads one to hope that in this situation the artifact might be removable. 

1.4. Objectives and scope. However, the explanation in the previous paragraph 
is informal in the sense that it ignores the singular behavior of the operators at 
the intersection An Cj, j = 0, 1, 2, and the intuition above needs to be justified. 
As an initial step, one should understand the mapping properties of operators in 
IP''' {A, Cj). The primary goal of this paper is to obtain sharp Sobolev estimates 
for operators in /^^'(A, Cj), j = 0, 1, 2. 

Although there exists a microlocal normal form for a folding canonical relation, 
C, obtained by applying suitable canonical transformations on the left and right 
[28] , when applied to A these same transformations will turn A into some unknown 
canonical graph. Thus, there is no such normal form for pairs such as (A, C) 



satisfying (i,ii,iii) in §S 1.1 For this reason, our results will be limited to the pairs 
(A,Cj), j = 0,1,2, described above. In other words, our results apply to general 
operators with the same cubic geometry underlying H,J7q; the geometry arising 
from F*F for any FIO with a folding canonical relation, including the single source 
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forward operator; and the geometry for the normal operator in the marine seismic 
imaging problem. In these cases, one has good control over the phase functions in 
the oscillatory representations of the operators. 



1.5. Main result. We show that, for the canonical relations that we consider, the 
following Sobolev estimates hold: 

Theorem 1.1. Let A e /P''(A,Co) or A e /P''(A,Ci) or A € /P''(A,C2), where 
the canonical relations Cj are as above. Then A : H^^^p — > ^loJ /"'^ 



(1.7) 



Remarks: 



p+i+e, / = -i,Ve>0, 
p -r ^ , 2 ^ ^ 2 ' 



1. Here H'^ , i?comp ^-nd iJf^^ denote the standard L^-based Sobolev space and its 
compactly supported and local variants. Thus our analysis of the Schwartz 
kernels can be restricted to compact sets in the spatial variables. 



2. Away from A n Ci, which is also the fold surface of Ci, the inclusion (1.6) 
holds, and /p+'(A\ Ci) and I^{Ci \ A) satisfy standard Sobolev space estimates 
associated with local canonical graphs. These estimates are always more regular 
than the ones mentioned in the statement of Thm. |1.1| Thus it suffices to only 
consider operators A supported microlocally close to An Ci. Keeping in mind 
both the and FIO natures of /P''(A,Ci) and the loss of g derivative for 

FIOs associated with folding canonical rel ation s [55], in general one certainly 



needs to take r > max(p + l,p + g). Thm. 1.1 can therefore be interpreted as 



saying that, when the strengths of the ^I^DO and FIO parts of A are sufficiently 

i<l<h 



close, specifically —l<l<h, there is a further loss due to their interaction. We 



will see in S|3] that this loss can in fact occur and Thm. |1.1| is in general sharp. 
For I < i close to | , this type of behavior is already present in the estimates of 
Greenblatt [TS] for fractional integrals along families of curves. 

For comparison, it is natural to investigate the analogue of Thm. |l.l| when C is a 
local canonical graph intersecting A cleanly in codimension one. We address this 



in Thm. 4.1 Singular Radon transforms over hypersurfaces satisfying the rota- 
tional curvature condition are of this general type but with higher codimension 
intersections. A particular case was first considered by Geller and Stein [T2| and 
the general class is due to Phong and Stein ; see also Cuccagna [7| , which al- 
lows for degenerate canonical relations. Different proofs of L^-boundedness were 
given in [TH1[37]. Our proof of Thm. |4.1| is essentially an adaptation of arguments 
in [18] using the methodology of parabolic cutoffs, cf. |27], whereas the more 
degenerate geometry of Thm. |1.1| requires a combination of a parabolic-type 
decomposition with a two-index dyadic decomposition inspired by the almost 
orthogonal decompositions of degenerate oscillatory integral operators due to 
Phong and Stein [551 
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1.6. Layout of the paper. In f|2j we recall the basic theory of Fourier integral 
operators associated with a single smooth canonical relation, as well as the paired- 
Lagrangian operators associated with two cleanly intersecting canonical relations. 
The oscillatory representation of the two classes of operators under consideration, in 
particular the normal operators for the linearized marine seismic imaging problem in 
the presence of fold caustics HT] is also reviewed. In f|3j we prove that the regularity 
exponents in Thm. |l.l| are sharp, in particular s how ing that the loss of derivatives 
for IP''' {A, Co) in the range — ^ < / < | in Thm |l.l| is unavoidable. As a warm-up 
to the proof of Thm jl.ll in f|4] we give an analogous result in the case when the 
canonical relation is a canonical graph; here, a standard parabolic decomposition 
suffices to establish the Sobolev boundedness. In f|5]and ^we introduce, for the 
classes /^''(A,Co) and /P''(A,C2) respectively, the phase space decomposition of 
mixed type which is the main technical feature of the paper. The proof of Thm. 
|1.1| is given in S|7j Finally, in fjS] we examine the Hilbert transform H along the 
cubic in more detail. We show that, although elliptic in a naive sense, H is not 
microlocally invertible using operators bounded on L^, indicating the difficulty of 
constructing parametrices for these kinds of operators. 

We thank the referee for suggestions concerning exposition and references. 

2. Background 

In this section we develop the notation and terminology needed in the sequel, and 
in particular provide an explicit representation for the operators of interest. We 
refer the reader to [23l El for the theory of classical FIOs, [29l [22l [30l [E] for 
material on the /^'' classes and parabolic cutoffs, and 11, §3] and the references 
there for further details on the relevance of these techniques for linearized seismic 
imaging. However, for convenience we briefly recall the salient facts in the form 
needed here. 

2.1. Preliminaries. Let X and Y be smooth manifolds and T*X \ 0,T*Y \ 
their cotangent bundles (with the zero-sections deleted), equipped with the canon- 
ical symplectic forms ujt'y, ujt*y- If A C {T*X \ 0) x {T*Y \ 0) is a conic 
Lagrangian submanifold with respect to ijJT'{XxY) — ''^l'-^T'X + '^*b}^t*y, then 
A' : (x, ^; y, — 77) G A} is called a canonical relation. For example, 

the conormal bundle of any smooth submanifold Z <zY y.Y is a. canonical relation. 
In particular, iiX = Y, then A = iV* A^, C (r*y \ 0) x {T*Y \ 0) is the diagonal 
relation. 

If A is a Lagrangian, meMandO<5<l-p<p<l, then /"^(A x r,A) 
denotes the space of Fourier integral distributions of order m and type (p, 5) whose 
wavefront set is contained in A (see ^ p. 97]), while if C is a canonical relation, the 
associated class /™^(A, Y; C) of FIOs of order m and type [p, 5) mapping £'{Y) — > 
2?' (A) consists of those operators A whose Schwartz kernels are Fourier integral 
distributions Ka € /"^(A x Y; C) This is abbreviated to I^giC) if X and 

Y are understood, and if p = 1 and 6 = 0, the subscripts p and 6 are omitted. If 
ttljTTr '■ C — > x*Y \ are the natural projections on the left and right, then a 
canonical relation C is a local canonical graph if and only if one of the projections 
(and hence the other) is a local diffeomorphism; of course, this can only happen if 
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dim X ^ dim Y. In this case, A e /"(C) =^ A : H^ompiY) — > H^oc"'^^)- 
In particular, if C = A, then /™^(A) is the class of pseudodifferential operators of 
order m and type {p,S), ^'"^(i^)- 

If C fails to be a canonical graph, then the critical points of the two projections 
are the same, and dim Ker [Dttl) = dim Ker {Dttb) at all points. One says 
that C is a folding canonical relation |28| , also called a two-sided fold |16| , if both 
projections have only Whitney fold singularities; in this equidimensional setting, 
that is the same as ^i^o singularity in the Thom-Boardman notation [13^ . A well 
known result of Melrose and Taylor [28] states that, for such a C, there is a loss of 
H derivative: 

6 

(2.1) A e /"(C) A : Hl^^p //^-^-J, Vs G M. 

Remark. Note that, combined with the composition result of |32J[TU], this gives 
an example exhibiting the loss in Thm. |l.l| for / = 0. In fact, let C be the folding 
canonical relation in the single source fold -caustic seismic geometry [35] , and let 
A G /™(C) be properly supported. By (O) the operator A* A G /^'"'''(A, Ci), with 
(A, Ci) as in Thm. A. Since the loss of ^ derivative for A is in general sharp, so is 
A* A : — !■ //''^2™~5, already showing the necessity of the loss of derivative 
for / = 0. 

2.2. Paired Lagrangian spaces. For an exposition on clean intersection theory 
and Fourier integrals, we refer the reader to [TS] [9]. Classes of distributions as- 
sociated with two cleanly intersecting Lagrangian manifolds were introduced by 
Melrose and Uhlmann [55] and Guillemin and Uhlmann [52] (see also [TH] for fur- 
ther relevant discussion). For our purposes, the definition of this class can be given 
in terms of multiphase functions 'HUI and symbol- valued symbols [551 HH] : which we 
now describe. 

Definition 2.1. Given a manifold X and a positive integer M > 1, let T be a cone 
inX X M*^ \ {0}. 

(i) A function ip = ip(x;0) G C°°(r) is a phase function if it is homogeneous of 

degree 1 in 9 and has no critical points as a function of {x, 6). 
(ii) A phase function is said to he non-degenerate in F if 

dgip{x] 6) — 0, {x, (?) G F the collection of vectors 

— — ■ j — ' ' ' , M j is linearly independent. 

(Hi) Let K he a conic C°° suhmanifold ofT*X\0. A nondegenerate phase function 
tp is said to parametrize A if 

A = {{x, d,^{x- 0)) : defix, 9) = 0, {x, 9) G F}. 

Definition 2.2. Let (Aq, Ai) be a pair of Lagrangians in T* X \ that intersect 
cleanly in codimension k. Let Aq G Aq H Ai and F C A x (M^ \ 0) x M*^ an open 
conic set. A multiphase function <j) parametrizing the pair (Ag, Ai) is a function 
(t)[x\9\a) G C°°(F) such that 

(i) 4'q{x] 9) :— (j){x; 9; 0) is a nondegenerate phase function parametrizing Aq in a 
conic neighborhood of Xq, and 
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(ii) (f)i{x; [6, a)) :— <p{x; 9; a) is a nondegenerate phase function parametrizing Ai 
in a conic neighborhood of Xq . 

Example. It is known [291I22| that any two pairs of cleanly intersecting Lagrangians 
are microlocally equivalent. One can thus consider the model pair (Ao,Ai) in 
T*M" where Aq and Ai are the conormal bundles of {x — {xi, ■ ■ ■ ,x„) — 0} and 
{xk+i = ■ ■ ■ = Xn = 0}, resp., so that 

Ao-To*M"-{(a;,0 :a; = 0}, 

Ai = {{x, £,) : Xk+i = ■ ■ ■ = Xn = 0,^1 = ■ ■ ■ ,(,k = 0}. 

Then Lp{x;d',a) = x ■ {a, 9'), with 9' e R"-'' \0, cr E M.'' , is an example of a mul- 
tiphase function parametrizing (Ao,Ai). In this paper, we will only be concerned 
with the case fc = 1, i.e., codimension one intersections. 

Definition 2.3. The space S'^'^ (X x (M^ \ 0) x M) of symbol-valued symbols of 
orders p, I, is the set of functions a{x; 9; a) G C°°(A" x (IR-'^ \ 0) x K) such that, for 
every relatively compact K (~ X , non-negative multi-indices a G 1^ , /3 G Z and 
7 e Z , the following differential estimates hold: 

(2.2) \d^d^d2a{x;9;a)\ < C„,^,^,K(e, (a)'"'', 

for all {x,y) G K, where (^) = (1 + |^P)^, etc., throughout the paper. 

The quantities sup^^g^„{9,ay°'\-P{a)f^-^\d^d^d2a{x;9, (t)| are referred to as the 

seminorms for the class S^'K If 16*1 > \a\ on the support of a, then we say that 
the phase variable 9 is dominant. 

We next define the classes of generalized Fourier integral distributions associated 
with an intersecting pair of Lagrangians. 

Definition 2.4. // Aq, Ai C T*{X xY)\0 are smooth, conic Lagrangians inter- 
secting cleanly in codimension one, then the space of generalized (or paired La- 
grangianj Fourier integral distributions of order p, I g M associated to (Ao,Ai), 
denoted /^''(X x y;Ao,Ai), is the set of all locally finite sums of elements of 
/P+'(Ao) + /P(Ai) and distributions of the form 

(2.3) u{x,y)^ J e"^^'-''y''^'''''>a{x,y;9;a)dcrd9, 
where a £ S'P''(X x F x (M^ \ 0) x M), with 

, , ~ N +1 dim X + dim Y , ~ 1 

(2.4) P = P + l + ^ ^ , ' = -^"2' 

and (j){x,y;9;a) is a multiphase function parametrizing (Aq, Ai) on a conic neigh- 
borhood of a point Aq G Aq H Ai . 

Remark. Aside from an adjustment in the orders, the paired Lagrangian spaces 
are symmetric 01, 22\ in Aq, Ai: 

(2.5) F-' (Ao, Ai) = (Ai, Ao) D P+\Ao) + F(Ai). 

Finally, we define the classes of generalized (or paired Lagrangian) Fourier integral 
operators which are the subject of this paper. 
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Definition 2.5. (i) //Co,Ci C {T*X\0) x (r*y\0) are smooth, conic canonical 
relations intersecting cleanly, then /^^'(X, Cg, Ci) denotes the set of operators 
A : £'{Y) — > V'{X) with Schwartz kernels KA{x,y) £ IP^\X x Y-C'q,C[). 

(ii) In particular, if C G {T*Y\0) x {T*Y\0 is a canonical relation intersecting the 
diagonal relation A cleanly, then the members of the associated class /^''(A,C) :— 
jp,i-(Y, Y; A, C) are referred to as pseudodifFerential operators with singular symbols 



The spaces X and Y are suppressed if clear from the context. Equipped with the 
definitions above, we are now ready to describe the operators of interest. 



2.3. Oscillatory integral representation of operators in IP'\A,Ci). The op- 
erators we consider will be of the following two forms. 



2.3.1. The cubic model. For the cubic model, the operator A lies in the class 
jp,i(Y^ y. (^^-j^ where Y is an open bounded subset of M^, and Cg is the conormal 
bundle of the cubic {x2 — 2/2 = {xi — yi)^}- The kernel of A e /^''(A, Co) can be 
written as 

(2.6) KA{x,y)= f ^'^'>^^^y^'^a{x,y;02;ei)d9, 

where the symbol-valued symbol a G S^^^''-^^ {Y x F x {Rg^ \ 0) x M^J, and the 
multi-phase function 0o is given by 

(2.7) Mx, y, 0) = {xi - yi)ei + {x2 -y2~ {xi - yif)e2. 

The reader may verify that 4>a{x, y; 62', 0i) is a multi-phase function parametrizing, 
as in Def. 2.2 the cleanly intersecting pair of canonical relations (Cq, A), in that 
order. The dominant phase variable is ^2- 



2.3.2. Normal operators for seismic imaging. We now turn to the operator classes 
jp,i(Y^ y. j _ 2^ 2, where (A, Cj) are the geometries that arise in the analysis 

of the normal operators F* F in seismic imaging, as described in the introduction. 
Here F is a bounded open subset of M", n > 3 (n = 3 for the physical problem). 

For the single source problem, the Melrose-Taylor normal form [3H] for folding 
canonical relations, such as the one to which F is associated, gives rise to a partic- 
ularly simple normal form for the pair (A,Ci), namely (A,Co x A7^*j^n-2). Th is 



allows for the estimates for /^^'(A, Ci) to be proven as for /^''(A, Cq); see §^ 7.1 



In contrast, the marine data geometry is complicated by the fact that there is 
only an approximate normal form for folded cross caps. It was shown in [11| that, 
microlocally near any point in AnC2 the pair (A, C2) can be parametrized, in that 
order, by a multi-phase function of the form 

(2.8) </.(x, y; e; p) = (x - y) . e -f ^ " i^!^±^a - p) , on 

(2.9) {\^i\>lm + \p\)}- 
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Here, P — P{xn,y,^') is an unknown smooth function on M x 
homogeneous of degree one in ^' = (^i, • • • , ^„_i) and satisfying 



r X (] 



1)71 — 1 



\0), 



(2.10) 

and 

(2.11) 



PI 



{x„=y„} 



0, VF|{, 







p\\ « icll 



microlocally, i.e., one can arrange for the left hand side of (2.11) to vanish at any 



chosen base point, and thus be small nearby. Here, is a dominant phase variable. 
Note that we have translated the formulas in [11], which were written in n — 1 
dimensions, to an n-dimensional setting. 

However, in order to handle this geometry in a way similar to the cubic model 
(A, Co), we want a multiphase function that parametrizes the pair (A,C2) in the 
reverse order. A prescription from |30j . which we now describe, allows us to compute 
such a phase function explicitly. We observe that A' is the flowout from C'2 D A' by 
the Hamiltonian vector field Hq on T* (M" x R"), where q{x, ^, y, rj) — (x„ — ?/„)^i. 
Starting with ^Q(x,y;^, p) = </) as in (2.8), which parametrizes C2, we solve for 
'4'{x,y;^,p;s) satisfying 

dip , .d^ ~ 

^ = (x„-y„)^, V^U.o = Vo, 

and find that 

?1 4 4i 

Finally, letting s = this yields the multi-phase function 

(2.12) V'(a;,y; (^,p);cr) = (a; - y) • C + - ^{xn + Vnfp- + {^n ~ yn)cr, 

4l 4 4i 



which par ame trizes (C2 , A) in that order according to Def. [2^2] ByDef. 2.4 and the 
relations and the kernel of an operator A e /P''(A, C2) /p^'(C2, A) 



has a representation, modulo /p+'(A) + I'^{C2) of the form 



a{x,y;^,p;(7)d^dpda, 



(2.13) KA{x,y)^ J 

with ae SP-^2^^-i (y xY X (M^'+i \ 0) x M^), and supp(a) C {|cr| < 



2.4. A technical tool: the method of stationary phase. In the subsequent 
sections, especially SjS] and { 6] we will need to estimate repeatedly the oscillatory 
integrals that arise from ( 2.6 ) and (2.13) as the Schwartz kernels of the compositions 
B*A and BA* , where A,B e /''■'(A, Cj). We will use the method of stationary 
phase to describe the asymptotic behavior of integrals of the form 



/(a. A) 



a e 



as A — > 00; for the sake of completeness, we state it in the form needed. The phase 
function ip is assumed to be real-valued and smooth, while the amplitude function 
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g is assumed to be smooth of compact support in the variables {w,a) e M" x 
with the following growth restriction: there exists some < \ such that 

(2.14) g = 0(A™+''l"l) asA^oo, uniformlyin(w,a). 



We state here the principle of stationary phase in a form that we need. 

Lemma 2.6 ([51, p. 14). Let (p{w, a) — ^{Q{a)w, w), where Q is a real n x n non- 
singular, symmetric matrix depending continuously on a, and let g be an amplitude 
satisfying {2.14) above. Then I [a, X) has an asymptotic expansion 

K^a, A) ^ (^) " |det(Q(a))r*e^^^"(«('')) £ iEMl^X'^ 



as X ^ oo, uniformly in a. Here 



which is a second- order partial differential operator in w. 



In the sequel, justifying an application of stationary phase will involve identifying 
an oscillatory integral with (quadratic) phase Xip and amplitude g, and verifying 
that (a) the Hessian of Xip evaluated at its critical point is nonsingular, and (b) 
the quantity (X~^RYg evaluated at this critical point decays exponentially in r. 
Also, while the notation ~ was used in the above lemma to denote an asymptotic 
expansion of /(a, A), in subsequent applications we will use to represent just the 
leading order term of the asymptotic expansion, with the understanding that the 
latter also provides the sharp size estimate for I{a, A) for A ^ 1. 



3. Sharpness of Theorem 11.11 



We start by showing that, at least for /P''(A, Co), the estimates in Thm. [13] cannot 
be improved in general. For I < the inclusion TP{Cq) C /P''(A,Co), cf. ( |2.5[ ), 
implies that the Sobolev mapping index r of /^''(A, Cg) must be at least as large 
as that of (2.1 1 for /^(Co), i.e., r > p + g. On the other hand, the fact hat 
^fP+'(]R2) = /pT'(A) c /P''(A, Co) implies that we must have r > p + L 

In the critical interval — | < ^ < ^ , we demonstrate the optimality of the estimates 
of Thm. 1.1 for the fractional integral operators along {t,t^ 
define, as in ( 1.3 ), 

(3.1) 



For 



< I < 



Ji+U{x)^f*[5{x2^x\) 



X{xi)dxi 



\Xl 



where x e C^{m.), x{t) 



1 for \t\ < 



Writing the kernel of the operator 



as an oscillatory integral, it follows from the discussion in §^ 
/"3''(A,Co). Thm. [li| therefore implies that Jl_^_l : 
for r 



21-1 
6 



2.3.1 that Ji^i e 

That this cannot be improved is seen as follows. Since the Fourier 



multiplier of J'l 



1+ 



1 IS 



(3.2) 



x(t) dt 
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it suffices to show that 

sup \m{OM\^ > c> 0. 
Using the substitution u — ^2^'^, du = {3^l^'^u'^^^)^^dt, one has 



(3.3) 



, 21-1 



Xjt) dt 



\U\ 6 



Since the improper integral /e ^ du converges to a nonzero value, the 



same is true for the integral in ( |3.3p . 

For ^ = — i, we replace |a;i|~° by log|xi| in (3.11, which still defines an operator 



j7o G / 2> 2 (A, Co) since log is conormal of order —1 for {xi = 0} on M. The 
corresponding logarithmic divergence of |m_i (0,^2)1 • 1^21^ then implies that one 
cannot in general eliminate the loss of e derivatives when I = — \ in Thm. 



1.1 



4. The canonical graph case 

As a warmup for the proof of Thm. |1.H and to provide a background for the de- 
compositions that are needed, we first prove an analogous result under the optimal, 
nondegenerate microlocal geometry, namely when C is a canonical graph. While 
the intersection C H A can be quite singular, we restrict ourselves to the config- 
uration closest to that of Thm. |1.1| namely when C intersects A cleanly, say in 



codimension k. For example, if 7 in ( 1.1 1 is replaced with the parabola (t, t ), then 



that % e 7^2 '2 (A, C) for a local canonical graph C, although for families of curves 
in higher dimensions, the canonical graph condition must always be violated some- 
where. We mention in passing that singular Radon transforms associated with a 
family of fc-dimensional surfaces in M" belong to the class /~2'2(A,C), where C 
is a local canonical graph under the rotational curvature condition 36J, and with 
A intersecting C cleanly in codimension k |18j . Estimates for such operators are in 
[5n [12 [Ml [m EH ES] ; see also [7] for worst-case estimates in degenerate cases. 

The following result is similar to Thm. jl.ll but much easier to prove due to the 
nondegeneracy of C . For C a conormal bundle, related estimates were found by 
Seeger and Wainger P3] . 

Theorem 4.1. Let Y he a manifold of dimension n. Let C C (T*y\0) x (T*F\0) 
he a local canonical graph intersecting A cleanly in codimension fc, 1 < fc < n — 1. 
// A e /f -'(A, C), then A : H^,,„,^ ^ Ht^J for 



(4.1) 



p-he, ?==-|,Ve>0, 
p+^, -|<?<|, 
[p + l, l>l 



Proof. We prove the theorem using a standard parabolic cutoff. Parabolic cutoffs 
in various guises have been used often in microlocal analysis; the argument given 
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here is a direct adaptation of that in |18 , bu t related ideas are also in [311 [371 127], 
Since A e /P''(A, C) = IP+''-'{C, A) by we can assume that A has kernel 



(4.2) 



KA{x,y) 



where is a multiphase function parametrizing the cleanly intersecting pair (C, A) 

and a & SP-^'^-iiY xY X (M^ \ 0) x M). Decompose 



a = X 



a+\l-x 



a := flo + ai, 



with A = ^0 + ^1 the corresponding decomposition of A. Since (cr) > c{9)^ 
on supp{ai), the differential estimates on ai = ai{x,y; (9, a)) ensure that it is a 
standard (non product-type) symbol, 

(4.3) 



,0 



(Y xY X (]R^+'= \ 0)) 

(y X y X (M^+'= \ 0)) , 

Since <f>{x,y; {9,(t)) parametrizes A, this implies that 



«1 € S p+;_iV^^ 



I > 



(4.4) 



^1 e *i 1 



max(p+2^,p+0 



We now turn to Aq, whose Schwartz kernel is of the form (4.2), but with the am- 



plitude a replaced by oq. Recalling that (j){x,y;d,0) parametrizes C and reasoning 
as in [m Prop. 2.1], let us integrate out a and write 



KAQ{x,y) 



e''*'^'='y'^'°^bo{x,y;e)d0, where 



boix,y;e) = / e*(*(-'^^'^'-)-*(^'^^'^^°))ao(x,y;0,a)da. 



Keeping in mind the symbol estimates for ao and the homogeneity of (j) in {0,a), 
the argument in |18j carries over with minor modifications to yield 



(4.5) 



bo & Si 1 



and hence Aq e Ii ^ 



-(C). 



Here, and throughout, for t £ M, we set 
(4.6) 



t ift>0, 
if t < 0, and 

any e > for i = 0. 



Combining (4.5), (4.4) it follows from the Calderon-Vaillancourt theorem and its 
extension to FIOs of type (5,5) associated with local canonical graphs [51[TH] that 



A : m 



for 



, 1,, k. 2l + k 
r > max [p + - (1+ -^h^P^ ^,P + l 
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which yields (4.1 ). The calculations of sjsjare easily modified for the curve {t, t^) 



to 

show that Theorem |4.1| is sharp, at least for fc = 1. □ 

When C is a folding canonical relation, more involved decompositions of A, re- 
quiring additional knowledge of the multiphase functions, are needed to prove the 
analogous estimates. and f|7] are devoted to dealing with the cases when C 

has one of the forms described in §92.31 



5. Mixed Parabolic-Phong-Stein Decompositions for /^'''(A,Co) 

In the previous section, we described and applied the method of parabolic cut- 
offs in the context of FIOs associated with canonical graphs. There is also a 
method introduced by Phong and Stein [38t l39] for proving estimates using 
dyadic decompositions in phase space; this, along with variations, has proven to be 
a powerful technique for controlling oscillatory integral operators and FIOs with 
degenerate canonical relations [HI [51 [SI [T7] . In the present context, our analysis 
of operators in /^'''(A,Cj), j = 0,2, will combine these two decomposition strate- 
gies. In this section, we describe this for the class /''■'(A, Co), where we recall that 
Co = N* {x2 — y2 — {xi — yifY- Appropriate modifications of these decompositions 
for /P''(A, C2) wiU be dealt with in s|6| 



5.1. Components of the operator. Let {xj : j > 0} be a non-homogeneous 
dyadic partition of unity on M such that 



(5.1) 



supp (xo) C {|t| < 2}, 

supp (x,) C {2^-1 < 1^1 < 2^+1}, J > 1, and 
|X-"^| < C„2-"^ foranm>l. 



Setting Xjk{0) = Xk{di)Xj{(^2), we note that {x^fcjO < k < j < 00} is a bounded 
family in S'^''^ with respect to the seminorms of (2.2). Furthermore, we set 

Xoo 1 ^ Xjkj 
j>Q k=0 

SO that Xoo S and is supported in {j^ij > ||^2|}- Given any A g /P''(A,Co) 
whose Schwartz kernel is of the form (2.6), its amplitude a(x, y; 62', 0i) decomposes 
as 



a{x,y]e2,0i) ^ \^^X]k{0) + XooiO)\ • a(x, y; 6*2; 6*1) 
\j>o fc=o / 
00 j 

:= ^^ajk{x,y;92;9i) + aoo(a;, y; 6*2; 6*1). 

j=0 k=Q 

Thus {ajk} is a bounded family in S'^'+i''"^ , with 

supp(a,fc) C {{x,y;e) : 2^-^ < |02| < 2^+\ 2'^-^ < < 2'=+i}, 
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and Coo G 'S'l o' • The decomposition of the amphtude a induces a decomposition of 
the operator 

oo 

j=0 0<k<j 

where the Aji^ and A^o are operators whose Schwartz kernels and muhi-phase func- 
tions are of the form described in §j ]2.3.1 (see (2.6) and (2.7)), but with the am- 
phtudes replaced by ajk and Uoo, resp. Thus, A^o <E /p+'(At.r2) = v1/p+'(M2-)^ 
{Ajk ■ < k < j < oo} C IP'' {A, Co). It follows from standard ^BO estimates that 



HZ 



s-ip+l) 



which is at least as smoothing as what is being claimed 
for the entire operator A in Theorem We therefore disregard Aoo in the sequel. 



For a fixed constant S e (0, 1) to be set later (we will eventually take (5 € [|, |)), 
we group the summands in A — A^o as follows, 

oo 

(5.2) A- A^ ^ Ao + ^Aj, where 

oo [Sj] j 

(5.3) Ao^Y.T.^i'^' E ^Jfc'J>l- 

j=0 k=0 k=[Sj] + l 

In this section we identify Ag as a standard FIO associated with the folding canoni- 
cal relation Co and record several almost orthogonality properties involving the Aj , 
Ajoo and Ajk- These facts will be used heavily in the proof of L^-Sobolev estimates. 

5.2. Bounds for Aq. 



2 ' 

where m — p + d{l + 5)41, and hence maps H^gj^p{M.'^) boundedly into Hi^^° (M?) 



Lemma 5.1. For any < 5 < |, the operator Aq lies in the class I^g g{Co), 
whe 
where 

(5.4) ro>p+l+5{l + ^)^. 



Here tr is as in {4-6)- 



Proof. The kernel of ylo is of the form (2.6), but with the amplitude replaced by 



00 [Sj] 
j=Q k=0 



4.1 



we 



We observe that supp (ao) C {(6*1) < 0(02)^}, and that ao satisfies the same differ- 
ential inequalities as x ((^i) (6*2)"*) •a(a;, y; ^2; 6*1), for some X S C'^(M). Integrating 
out 61 in the oscillatory representation of Kj^^ as in the proof of Theorem 
see that 

i^A„(x,y)= J e'''^---y^-(-'''~y'y"^ao{x,y;92)d02, where 
do{x,y;02)= I e''^^'-y'^''ao{x,y;02,0i)dei. 
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Figure 1. Decomposition of (6*2,6*1) space. 



It is easily verified from the support and dif fere ntiability properties of oq, using 
arguments similar to those used in Theorem 4.1 that ciq £ S-^^g ' , where m is as 



in the statement of the lemma. In other words Aq is an FIO associated with the 



two-sided folding canonical relation Co, of order m + 



— m and with 



an amplitude of type Since useful FIO classes need to be invariant under 

changes of variables and phase functions, we may as well c onsider (cf. [TH]) the 
type as being {1 — 6,6). Now, 1 



since 6 < -k, so (5.4) follows from 



□ 



5.3. Almost orthogonality. The next two lemmas deal with the almost orthog- 
onality among the components {Ajf^}. Both for clarity and for possible future use, 
we will treat the A* in A* A as the adjoint of a general operator B £ 1^ (A, Co). 

Lemma 5.2 (Almost orthogonality in j). Fix < 6 < 1, and orders p,p' , I, I' £ M. 
Given A e /P''(A,Co) and B £ I'p (A, Co), we decompose A and B as described 
in %\5.1 using the same value of 5 for both. Let {Aj^ '■ Q < k < j} and {Bj^. : 



< < j} denote the components of A and B respectively (see equations {5.2) and 



AN FIO CALCULUS: SOBOLEV ESTIMATES 



17 



Then for every N > 1, there exists a constant Cn > such that for all 
> 1, \j — j'\ > i, Sj < k < j , Sj' < k' < i' , the following estimate holds: 



(5.5) 

In particular, if \ j — j'\ > 3, > I, then 



R* A 

\-^]'k'^jk\\l^2^L2 



< Cn2 



-AT I 



(5.6) 



\B*,A, 



< Cn2' 



-Ni 



U-j') 



The same estimates also hold for Bjik'A*^. and BjiA*. 

Proof. We only give the proof for B* A, that for BA* being similar. The Schwartz 
kernel of the composition B*n,,Ajk has the form 



(5.7) KB.,^,A,,{x,y) = j e'-^^-'y''-fi-fi'^b,,k'{z,x-9')ajk{z,y\e)dede'dz, 

where the phase function $ is given by 

(5.8) $(x, y, z- 9, 9') = (zi - yi)0i - (zi - xi)9\ 

+ {Z2 - 2/2 - [zi - yif)92 - {Z2 - X2 - {zi - Xi)^)9'2. 

We note that dz2^ ^ 92 — 9'2, so that for \j — j'\ > 3 and on supp (bj'k/ ■ CLjk), 

\dz $1 > {2™'^^'^^'^'^~^ — 2™'"(^'-'')+-'-) > 12™^''(-''J') 

Integrating by parts N times in Z2, integrating in all variables, and applying the 
support and differentiability properties of Oj^ and hj'k' throughout we obtain 



(5.9) KB'.,^,A,Ax.y) 



2(P+5)j + (p'+5)/ + ('-3)'=+(''-3)'='-^"iax(j,i')^ 

Integrating this over the compact support in either x or y yields the estimate 



sup 



dy + sup 

y 



where M = 4+|p+i| + |y + i| + |Z— + — i|isa fixed constant. Schur's lemma 
then implies the same bound for ||i3*,j.,yljfc||i2_>.^2. Since N is arbitrary, the proof 
of (5.51 is complete. The inequality (5.6) follows from (5.5) simply by summing in 
k and k' . □ 

Lemma 5.3 (Almost orthogonality in k). Under the same hypotheses as Lemma 
|5.^[ the following conclusions hold. 



(a) There exists a constant C > .such that for all j > 1 and [Sj] + I < k < j , 

(5.10) \\Ajk\\L^^L^ <C2'^P+h)j+(i'h)k^ 

A similar .statement holds for Bjk. 
(h) Now suppose 8 > ^. Then, for any N >\, there exists a constant Cjy > such 
that for all fc, k' > 1, with \j - j'\ < 2, |fc - k'\ > 3, [Sj] + 1 <k < j and 
[Sj'] + 1 < k' < j' , the following estimate holds: 

(5.11) \\B*,k'Ajk\\L^^L^ < CAr2(P+P'+i)j'+('-3)fe+(''-3)'='2-l'^-'='l^. 
The same estimate holds for Bji^iA*^. 
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Proof. For any choice of indices j,j',k,k' with [Sj] + 1 < k < j and [6j'] + 1 < 
k' < j' , we make a preHminary simphfication of the integral kernel Kgt^^^j^^^, 
that will be useful in the proof of both ( [a]) a nd (|b|. Going back to the oscil- 
latory representation (5.7) with $ as in (5.8), we perform stationary phase in 
(62,22) on this oscillatory integral. This step is justified by Lemma 2.6 with 
w — (6*2,2^2), a — {zi,x,y,9,0'i) and Xtp = $ in the notation of that lemma, so 
that iX-^RY[bj,k'ajk] = {dg,^^J''[bj'k'a-jk] 0(2"^'^) < 1. We are thus led to the 
estimate 



(5.12) 



Kb- A, y) 



where z — {xi,X2 + {zi — Xi)^) and 

(5.13) $ = (zi - yi)9i - {zi - xi)9[ + {x2 - + {zi - xif - (zi - yif)92. 

Several integration by parts will now be performed on the integral in ( |5.12[ ). In 
the subsequent computation, > 1 will denote an arbitrarily large integer and 
Ujj'kk'{x,y, zi;92,9i,9[) any function satisfying similar size and differentiability 
estimates as bj>k' {z, x; 92, 9[) ■ ajk{z, y; 9), possibly with different implicit constants. 
The value of N or the exact functional form of Ujjikk' may vary from one occurrence 
to the next. 



For part ([a|, we set j — j' , k — k' , integrate by parts 2N times in (5.12) with 
respect to each of the variables 6*2 , 6*1 , 9[ , finally obtaining 



K 



B;,A,Ax,y)\<CN / [l + (2'=(zi-yi))'] " [l + {2''{zi-x,)) 



-N 



X [I + (2^X2 - y2 + {zi ~ xif - {zi - yif)f] ^Uj,kkdzid9[d0. 

Fixing X in the above expression, and integrating out the variables y2, yi, zi, 9[, 
in that order yields the estimate 



(5.14) 



sup 



ifB-A,.(^,y)U2/<C2(P+^''+i)^+('+''-i)* 



The same estimate is obtained if the roles of x and y are interchanged in (5.14). 
Schur's lemma then gives the following estimate for the operator norm: 



\B*kAjk\\L^ 



»L2 



< 



sup 



2 






sup / 




. y J 



\K 



{x,y)\dx 



(5.15) < C2^i+i'-^)k+(.P+p'+i)o^ 

Substituting B = A, p = p' , I ^ I' yields the conclusion of part ([a]) of the lemma. 
We now turn to part ([b|, where we deal with the "off-diagonal terms" B*,j^,Ajk 



with \i — / I < 3 and \k — k'\ > 3. Without loss of generality, we may assume that 
fc > A;' + 3. Observing that 



{x^ ~ yi)e'* = [(zi - 2/1) + (xi - yi)]e** = (dg, + dg,)e 
an integration by parts in (5.12) in the 9i, 9[ variables gives 
[l + {2'^'{x,-y,))Y\KB',^,A,,{x,y)\<CNX 
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^t^(x,y,z-e,e[)^^^^^^^ (a;, y, z; 6, 9') dzi d9[dd 

On the other hand, 

4,l> = 9i-e[+ W2 ((zi -xif- (zi - 
= 01 - 302(4 

from which we obtain L(e**) = (0i -6'i)e**, where L := jd^^ -S^adg^ +36'2rfe, . We 
observe that j^i — 9'i\ > ^2'^, and that the loss from applying L* to the amplitude 
(6*1 - 9[)-'^Ujj,kk' is < C(l + 2J-2fe _^ < C'max(l,2J-2fe'). More precisely, 

for any m > 1, 



min 1,2 



9i - 6*;)™ 



'^jj'kk' 

= (01 -ei)"- 



Thus an integration by parts argument using the differential operator yields 
(5.16) [l + (2'='(xi-2/i))2]'^|AB;,^,A,,(a;,2/)| < C^x 
l+(min(l,22'='-^)(0i-0;))' 



-"jj'fcfc' (a;, y, zi; 6*2, 6*1, ^i) dzi d9 

We observe that the argument so far does not require any special choice of 5. 

i 

2 



We now consider two cases. First suppose that k' > |, so that min(l,2^'^ ^) — 1. 
Since j^i — 05.1 > > |:2''^ on the support of Ujjikk', the estimate in (5.16) implies 



\Kb', ,A,,{x,y)\dy+ / ,A,,{x,y)\ dx < 0^2^^-'^'' < Cn2~^^ 

for arbitrarily large N . Recalling that |fc — fc'| < j, the inequality above combined 
with Schur's lemma yields a stronger estimate than the one claimed in (5.11 ). 

Next we as sume that 2k' < j. Here we will integrate by parts twice more in the 
integral in (5.161, using the differential operators = 1 + {^de^) , r — 1,2, and 
keeping in mind that 

deAe'^) = i{zi - yi)e'*, de.(e^*) = i{x2 - y2 + [zi - xif - {z, - yif)e'^. 
We observe that 

de^Ujjfkk') = 2^^Ujj,kk', while 



Ujj'kk' 



(l + (22'='-J(0i -0i))2)A^ 



Ujj'kk' 

(l + (22fc'-J(0i-0i))2)A" 



(l + (22fc'-J(0i-0'i))2)JV 

where the last inequality follows from the assumptions k > k' > (since ^ > |) 
and li— j'l < 3. Combining all the arguments above we finally arrive at the estimate 

[1 + (2''^'(xi - y,))Y\KB'^,^,A,dx,y)\ < C^2('-^)'=+(''-^)'='+(^'+^''+i)^- x 
"l + (22^'-^(0i - ^'i))']""^ [l + (2^-2'=' 1^1 - yi 
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1 + (2^|X2 - 2/2 + {zi - x^f - (zi - y^f\y 



-N 



dzi dO[ d6 



Choosing N large enough and observing that 16*1 — 9'^\ > ^2 , we obtain 
\KB':,^,A,Ax,y)\dx 

N 

V 21 ~2 



/I 



(5.17) 



[l+(2^--^'='|zi-2;i|)^]-^ 



[l + (2^ (xi-yOf 



N 



dx2dxidzi dOidO'i d02 
1 + (2J>2 ~y2 + {zi - xi)3 - (zi - yi)3|)' 
< (72('"3)''+(''"3)'''+(''+p'+i)-''2~^^^'''+''~-'\ 

where the last step follows from the fact that the {x2, xi, zi,9i,9'i,92) integrals 
when computed in that order yield 2"-' , 2^*^ , 2^*^ , 2^^'^'' , 2*^ and 2-' respectively. 
By symmetry, the same estimate also holds for the integral with respect to dy of 



\t^B*,^,Ajkix,y)\, with X fixed. In view of the last line in (5.17) and the assumption 
that 5 > |, (5.11) then follows from Schur's lemma and the inequality 



2k' + k-j^ 3k' -j + {k - k') > {3k' - j') + {j - j') + {k- k') 

>0-2 + (fc-fc') 

since k > k' > Sj' > j' /3 and |j — j'\ < 2, completing the proof of part (|b]). □ 

Corollary 5.4. Let A e /^''(A, Cq). Then for j > I and S > |, the operator Aj 
satisfies the property 



(5.18) 



11^. 



< 



9 [p+i +<?(;- 



ifl<\. 



Further, if p + I < for I > \, or if p + ^ + 6{l - ^ < for I < \, then 
^ = A — — is a hounded linear map on . 



Proof. The proof of ( 5.18 ) is a direct consequence of (|5.11|) in Lemma 5.3 with p ^ p' 



and Cotlar-Knapp-Stein almost orthogonality lemma |44j . The second statement 
uses (5.181, Lemma 5.2 and almost orthogonality again. The details are left to the 
reader. □ 



6. Normal operators for seismic imaging 

We now descr ibe th e deco mposition of an operator A e /^"'(A, C2), where (A, C2) 

As in §5jlet us fix 5 £ [7, 1), and a dyadic partition of 



is as in 1.3 



and 2.3.2 



unity {xj '■ j > 0} on M satisfying ([5.1 1). Setting 

j 

j>0 k=0 
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we arrive at a partition of unity on ^K^'^^ \ {0} 

oo j c 



namely 



j=0 k=0 



Xo 



E 

3=0 



[<5il 

E+ E 

fc=0 k = [5j] + l 



Xjk + Xo 



Here Xoo G is supported on {{£,,p) < 2|cr|}, while {xjk : j > 0, < fc < j} is a 
bounded family (with respect to the seminorms) in S*"^*^, with 

supp(x,fc) C {2^-1 < {^,p) < V+\ 2^-1 < \a\ < 2'=+!}, so that 



(6.1) 



supp(^^X,fe) C {\a\ < C{C,p)'}, and 

j>0 k=0 

j 



j>0 fc=[<5i] + l 



Letting ajkix,y]^, p;a) = Xjkit P,(7)a{x,y;^, p;a) and = XooO gives rise to a 
decomposition of the amplitude a, which in turn induces a decomposition of the 
operator A. More precisely, 



A^Ao + Y, E ^^'^ 

3=0 k=[5j] + l 



A,, 



where Ajk, Aq and A^o are operators whose Schwartz kernels and multiphase func- 
tions are of the form given in §^2.3.2 (see (2.12) and (2.13l), but whose amplitudes 
are given by a^k, J2j>o J2o<k<[Sj] °-3k and aoo respectively. Then A^o G and so 
maps H'^ boundedly into if*"'' for r > p + l. It turns out that the components Aq 
and Ajk in this situation satisfy estimates and almost orthogonality properties 
analogous to their /P''(A,Co) counterparts with similar but more involved proofs. 
We record these facts below with the appropriate modifications. 

Lemma 6.1 (Bounds for Aq). For any < S < ^, the operator Aq £ si^^) , 
and hence maps iJ^Q,^p(M") boundedly into Hi^J"{R"), where m and rg are in 
Lemma \5.1[ 

Proof. The argument is identical to the one presented in Lemma |5.1| and involves 
integrating out a in the oscillatory representation of Kaq- The details are left to 
the reader. □ 

Lemma 6.2. Let A e /P''(A,C2), B e /p'-'' (A, C2), with decompositions {Ajk} 
and {Bjk} as above. Then the conclusions (5.5) and (5.6) oj Lemma 5.2 hold for 
the same set of indices j , j' , k, k' therein. 

Proof. The kernel of B*,f.,Ajk is 



(6.2) 



B'-,^,A,Ax.y) 



X bj'k' {z, x; ^, p; a)ajk{z, y; p; a)dzd£,dpdcrd(_dpda , 



where, from (2.121, 



* = {z - y) ■ £, - {z - X) ■ i + ^ - ^ + {Zn - yn)cr - {Zn - Xn)a 
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i {{Zn+ynfp-{Zn+Xnfp)-P{Zn,y,£,')f^+P{Zn,X,^')f-^, 



and — (^1, • • • , Cn-i)- In view of ( |2.9| ) and (6.1), the amplitude satisfies the 
support condition 

^ n(e,P>l<2||i|, 2^-'-i<|(|,p)|<2^-'+\ 2^'-i < |a| < 2^'+M 

I l(e,p)l<2|6l, 2^--i<|(e,p)|<2^+\ 2'^-i<|a|<2^+ij 
and the differentiability estimates 



(6.3) 



with C 



independent of j and k. 



Simply integrating the zeroth order bounds in (6.3) yields the basic estimate 



(6.4) 



j'k' ' 



Now, if i' < j — 3, then |fi — ^i| > c|fi| > c2^ on s\xp\i{bj'k'0.jk)- Noting that 
dz^^ = S,i — ^1, can integrate by parts N times in zi and then estimate as in 
(6.4) to obtain 



< C7Ar2-(^-*^)^' 

for some fixed M — M{p,l,n) > 0, since k' < j' < j and k < j. By choosing 
N sufficiently large compared to M, integrating the kernel in x and taking the 
supremum in y, or vice versa, we get the same upper bound, which establishes the 
desired norm by Schur's lemma. The case of j' > j ' + 3 is identical. □ 

Lemma 6.3 (Almost orthogonality in k). Let A, B, {Ajk} and {Bj/t} be as in 



Lemma 6.2 If the indices j, f , k, k' satisfy the same hypotheses as in Lemma 5.3 
then the conclusions 15.10 ) and (5.11) hold. 



Proof. The proof of (5.10) and (5.11 ) in this setting requires an initial preparation 
of the Schwartz kernel of B*,f.,Ajk] we simplify (6.2) using four applications of 



stationary phase. The first three applications are with respect to the pairs (p, ^n), 
{p,^n) and (z',^'). Of these the first method of stationary phase may be justified 
using Lemma [2. 6[ by setting w ~ {p,£,n), so that in the notation of that lemma 











-26 
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d^ndp' 



and {X'^Rnbyk'ajk) = 0(2-^, 



giving rise to a valid asymptotic expansion whose first term involves a factor 
det((3)^2 = ^ 2^. The other stationary phases are handled similarly, their 
justification being left to the reader. By the same argument as before, the second 
stationary phase contributes a factor of 2^ (from the determinant of the Hessian) 
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to the leading order term of the asymptotic expansion, while the corresponding 
contribution from the third is only a constant. The end result is 

(6.5) KBp^,A,ki^^y) ^ J e^'^^c^dznd^'dada, 

where $3 and C3 are both functions of (x, y; z„; cr, a), with 
*3 = (x' - v') ■ ^' + [zn - y„)cr - (z„ - x„)ct 

- ^ - Xn){Zn + Xn)"^ - {Zn " Vn){Zn + Vn)^) 6 

+ {Zn- yn)P{Zn,y,0 ~ {Zn~ Xn)P{Zn,X,£!), and 



C3 = X bj>k'{x',Zn,x;^', -(z„ +x„)^^i + P(z„, a;, ^'); (a;„ - z„ 



X ajk{x',Zn,y;^', -(z^ +y„)^Cl + P{Zn,y,S.'); (Vn - Zn)^l,(T) 

The amplitude C3 satisfies the same size and differentiability estimates as 2^^^ hj'ya 



'jk 



with the eliminated variables absent (see (6.3)). 

To estimate the norm of each Aj^, we set B = A, j = j' , k — k' , p — p' , I — I' and 
apply one final method of stationary phase in {zn,a) (once again justified by the 
existence of the large parameter \a\ > ^ 1), obtaining 



(6.6) 



KA^^Ajx,y)^ J e^'^'c^deda. 



Here $4 and C4 are functions of (x, y; cr), with 



(6.7) 



$4 



{x' - y') ■ + {Xn - yn)<7 + -.{Xn " yn){Xn + ynf£.l 



+ {Xn - yn)P{Xn,y,C), 

and C4 satisfying the same estimates as C3, i.e., 

(6.8) \d^,d^d2^yC4\ < c„^^2(2p+ib+(2'-i)fe-l"b-l/3|fc. 

Simple integration in ^',(7 in the range |^'| < 2^+^,|(t| < 2*^+^ yields the basic 



estimate 



Ka' A,Jx,y) 



+2lk 



but one may substantially improve upon this by integrating by parts, using the 
differential operators 

dji$4 = Xi-yi + ^{Xn-yn){Xn+ynf + {Xn-yn)PiAXn,y,£,'), 

(6.9) de'*4 = x" -y" + ix^-y^)P^.{xn,y,0, 

d<T^4 Xji ym 

where x" = {x2, ■ . . ,x„_i). Since d^,P{xn,y, C) is homogeneous of degree 1 — |/3| 
in and |^'| > 2^-\ we find that 

d 



(6.10) 



C4 

{d^.<^4 
da 



C4 



id^,^4 



l<i<n — 1, m>l, 



while 



(C4) = 2 ''£4, 
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where £4 denotes a function satisfying the same estimates (6.8) as C4, possibly with 



different imphcit constants, and whose exact functional form may vary from one 
occurrence to the next. Integrating by parts N times with respect to each of the 



variables (^1, . . . ,^„_i,cr) in ( |6.6[ ), and combining (6.91 and (6.10) yields 

<Cn 



Ka- 



2(2p+l)j + (2i-l)fc ^^/^^ 



from which one obtains 



l + 2J'|de3'4|)^(l + 2'=|d<,$4 

d^'dadx 



N ■ 



(l + 2J'|d^-$4|)^(l + 2'=|da$4 



< (^2(2P+l)j + (2'-l)fe, 



l«'l<2- 
|ct|<2'=+i 



^1 dS^'da 



In the second step above and in view of (6.9), one has to integrate first in x' and 



then in x„, while the last step uses the size of the (^', cr)-support of a,fc. The 
same estimate holds for J \Ka* A k{x,y)\dy, completing the proof of (5.10) for 

/P''(A,C2). 



In order to prove the analogue of (5.11), it suffices to show that for any > 1, 
|j-/|<2, |A:~fc'|>3, 



.11) sup 



Kf 



',,,A^Jx,y) dy + sup / KB',^,A^Ax,y) dx 



< CAr2(P+^''+^^.''+('"5)'=+(''-5)'''2"l''~'''l 



N 



For this we return to the representation (6.5). As in the proof of Lemma 5.3 



we will subject the integral in (6.5) to a large of number of integration by parts 



using several differential operators. However due to the complicated structure of 



the multiphase function (6.6) (compared to (5.13)), the differential operators and 
hence the resulting integrations by parts are more involved. We now proceed to 
describe each of these steps systematically. Throughout this discussion, as in the 



proof of Lemma 5.3 we will assume that fc > fc' + 3, and denote by Ujj'i^ki any 
function (with possibly different functional forms) satisfying the same size and 
differentiability estimates as C3. The value of the large constant N may also vary 
from one occurrence to the next. 

Step 1. Since 

(6.12) do.$3 = z„ - ?;„, da^3 = Xn-Zn, (d^ + da-)$3 = x„ - ?/„, 

with 



dcr{ujj'kk') — ^y-jj'kk' and d^^Ujj'kk') — '^Ujjikk', 



integrating by parts in (6.5) a large number of times using 1 + [j('io- + dcr)] and 
1 + iidcrf gives 

N 



(6.13) 



1 + (2^= (x„ - y„)) KB*,^^,A,Ax,y) 



(a;,y, 2;„;cr, a) 
[l + (2'='(a;„-z„))2]^ 



dzn d^' da da 
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The above expression permits a localization in the spatial variables Xm Um Zn, which 
will be useful in the sequel. Fixing x € C'o°([~2, 2]), with % = 1 on [—1,1], and 
introducing the partition of unity 

1 ^ x(2'='-'''(^n - Vn)) + (1 - x)(2'='''=''(a:„ - y„)) 



= X(2 



[xn ~ Vn)) x(2'='-'='^(a;„ - z„)) + (1 - x)(2*'-'='^(x„ - z„)) 



k' — k e 



+ {l-x)(2'''-''''(x^~yn)), 



we obtain the decomposition 



where each K.i is an oscillatory integral with multi-phase $3 , and amplitude of the 



form (1 + {2'''{xn 




+ (2'='(x„-j/„))')"'^0., where 


61 = 


X(2'='-'='^K 


- y„))x(2'' "'^ '(x„ - Zn))ujj^kk', 


62 = 


x(2'='-'^''(x„ 


- yn))(l - x)(2'''"''''(a;„ - Zn))ujj,kk 


63 = 


(l-x)(2'='-'= 


('^n Un^^'^-^jj' kk' ■ 



It follows from (6.13) that for any e > 



\ICiix,y)\ < Cn2-^^ for i = 2,3 and aU TV > 1, 



which is a stronger statement than the one required in (6.11 ). We therefore restrict 
attention only to ICi in the sequel, in which Xn,yn and z„ are further restricted to 
satisfy 

(6.14) \xn - z„| < 2l-'='+'='^ \xn - yn\ < 2l-'='+'='^ hence |2„ - 2/„| < ^2"'='+'='^ 



Step 2. Our next integration by parts (in /Ci) will involve z„, cr, tr and will 
exploit the disparity in the sizes of cr and a. For this, we note that 

dz„^3 = ^ "5- + ^ [2(a;„ - 2;„)(a;„ - y„) - (a:;„ - y„)^] + P{zn,y,£,') - P{zn,x,^') 

+ {Zn - y7i)Pz„{Zn,y,C) - (Zn " Xn)Pz^{Zn, X, £,') 
= a - [2{Xn - Zn)iXn - yn) - {Xn ^ J/n)^] 

3 3 

+ 2{Zn - yn)'^Pz^z„{yn,y,£.') - -^{Zn " Xn)'^ Pz„z^{Xn, X, ^') 
+ [Zn - yn)^Q{Zn,y,£.') ~ {Zn - X^)^ Q{Zn, X , ^') , 

where the last step follows by expanding P(z„, y, and Pz^ {zn, y, £,') (resp. P(zn, x, C) 
and Pz^ {zn, X, ^')) in a Taylor series in the z„ variable about Zn = Vn (resp. Zn — Xn) 
using (2.10). Here Q is a smooth function that is homogeneous of degree one in 
In view of (6.12), we find that 

L($3) = cr - CT + (z„ - yn)^Q{zn,y,^') - {zn - a;„)^Q(z„, X, where 
L = 4„ - I [2d^{d„ + di) - (rf, + 



- \Pz^zSyn.V.i')dl+^-Pz„zAXn,X,Odl 



L1+L2, with Li :— dz„, and 
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Here — — ^Pz^z„{xmX,£,') is a smooth function in all its arguments 



and homogeneous of degree one in The support properties in (6.141 and the 
homogeneity of Q imply that 



(6.15) 



|L($3)| > kl - |a| - C2^>„ - 2;„|3 - C2^>„ - x„| 



> 2 



fe-i o*; +1 



C2^ 



-3fe+3fc£ 



C2^ 



-3fc'+3/c'(! 



> 02" 



where the last inequality follows from k > k'+3 and fc+3fc'(l— e)— j > /(I— e) — ^ > 
j(i - e) — 2 >• 1, choosing e < |. 

We will integrate by parts a large number of times in /Ci using L, which is justified 



in light of (6.15). In order to describe the action of the differential operator on the 



amplitude, let us denote by Qi any (generic) function satisfying the same size and 
differentiability estimates as Qi, and observe that L2 is independent of z„, so that 
for any m > 1, 



^2 



(6.16) 



61 



L($3)' 



2 2-fe'(2-r)c) 

2J V^^ — ^— since fc > + 3 and T'Gi^ 2^61. 
^(L($3))™+-' 



On the other hand, 



61 



L(L($3))'"(1 + (2'='(x„ - z„))2)^J (^($3))™ [1 + (2fc'(a;„ - z„)) 



61 

(L($3))"+i 



1 

(L($3))'"+i 



(Zn - yn)'^Q{Zn,y,C) " - X,-,)'^ Q{Zn, X, Q 

[l + (2'='(x„-z„))2]^ 

^ 61 22fe'(.T„-z„) 

(L($3))"[l + (2'='(x„-z„))2]^+i 

' (L($3))'"[l + (2'='(x„-z„))2]^ 

' 2^{Zn~yn)^ei-2J{z^-Xn)^Ql 

Qi 22fc'(a;„-z„) 
(L($3))'"[l + (2'='(a;„-^„))2f 



where Q shares the same smoothness and homogeneity properties as Q, and hence 
QiQ — 2^Qi. In view of (6.12) and the factors involving (z„ — ?/„), (a:„ — Zn) in the 
expression above, one can follow up an application of L\ by another integration by 
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parts using da and d^, obtaining 



(L($3))"(l + (2'='(x„-z„))2) 



N 



dzn d^' da da 



.17) 



/ 



l + (2fe'(x„-z„))2 



l + (2fc'(x„-z„))^ 




ei 



(L($3))-+iy 

d2„ d^' da da 



dzn d^' dada 



Combining (6.16) and (6.17), we observe that a J-fold application of integration by 

am of oscillate 
,r form, name] 

2;3./(j,fe',p)Q^ 



parts in ICi using L results in a finite sum of oscillatory integrals with multiphase 
$3, whose amplitudes are all of a similar form, namely 



(6.18) 



(L($3))^ [1 + (2'='(x„ - z,,)yf [1 + (2fe'(a;„ - y^))Y ' 



with fi > J. It is easy to check using the restrictions 3k' < j and k' < k that the 
exponents /3i(j, fc', ^) satisfy the inequality /3i(j, fc', /z) < {fj,— l)k + k' . An induction 
in J, which we ask the reader to verify, then shows that in general 



(6.19) 



Pjij,k',fi)<i^i-J)k + k'j, J>1. 



We work with a general term in this sum, which by a slight abuse of notation we 
continue to denote by /Ci: 



.20) 



/Ci(a;,2/)~ J e**^:=,i(x, y,z„;^',cr, a) dznd^'dada, 



where Si is as in (6.18). 

Step 3. Our last integration by parts will be in the variables tr, a. Specifically, 
we note that 

- ^' ~y' +\ [i^n - yn){Zn + Vnf - {Zn ~ Xn){Zn + X^f] 6*1 

= {zn - yn)Pc{zn,y,^') " {Zn ~ a;„)P^' (z„ , ^') 

= Tl{zn,y,C){zn - Vn)^ - {zn - a;„)^7^(z„ , x, ^') , 

where the last step follows from Taylor expansion of P^/ (z„, y, ^') (resp. P^' (z„ , x,^')) 
in the z„ variable about the point z„ = ?/„ (resp. Zn = Xn) using ( |2.10[ ). Here 
TZ = (TZi,--- ,TZn-i) is a smooth vector- valued function that is homogeneous of 
degree zero in In view of (6.12), this implies that the vector differential operator 
P = (Pi,-- - ,Pn-i) = d^, -nizn,y,C)dl -nizn,x,C)dl Satisfies 

(6.21) P($3) =x' -y' + ^ [{Zn - yn){Zn + y-nY - [z-n ^ Xn){Zn + XnY] 6 1 . 
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In order to describe the result of a large number of integration by parts in JCi using 
P, we observe that P($3) is independent of {^',a,(j), so (in view of (6.20) and 
(6.18)) one only needs to understand the effect of P* on 8i(i($3))^^. For every 
1 <i<n-l, 



P 



ei 



©1 



(L($3))^- 



As in step 2, we exploit (6.12) and use another integration by parts with respect 
to cr, a to replace the factors involving (z„ — ?/„) and (a:„ — z„) in the expression 
above, thus obtaining 



2-iei 



3.22) 



d^'dada^ J dCdadae'"^ 

[QiA^n,y,adi + QiA^n,x,adi] 

[7^,(z„,y,e')rf' +7^^(z„,a;,e')4] ( 



2-^91 , A 2-'='(3-'-)ei A 2-'='(3-'')ei 



^0 W^^)) 



d^' da da, 



since Q^^Qi — ©i and TZiQi = 8i. Recalling (6.15) and observing that 
fc'(3 - r) + fc(l + r) > fc'(3 - r) + kr > 3k' > j 



we conclude from (6.22) after an easy induction that an m-fold integration by parts 
using Pi yields 



(6.23) 



ei 



i($3)' 



dS,' da da 



1 



{Pr^sY 



r) 



©1 d^dada 



fi+r 



where 

(6.24) lm{j,k',r)<-jm. 

This concludes the steps that involve various integration by parts. 



Step 4. It remains to combine the results of the previous steps to prove (6.11). 
Combining ( 6.20 ) with ( 6.23 1 , choosing to = M ^ 1 and continuing to call a generic 
term in the resulting sum by /Ci, we find that 

2^j(i.'^',A') 



K,i{x,y) < Cn 



[l + {2>^'{xn-yn))'- 



'"^^QidzridS^'dada 



N 



(L($3))'^+'' 1 + (2fc'(z„ - x„))^ nr=i [1 + 2-7M0-,fc',.)(p^$3)A/] 
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To prove (6.11) we compute 



sup f \ICiix,y)\dy < C2'^v+p'+i)3+(i-h)k+{i'~\)k' 2pjU,k' ,t.)-kt.^ 

X J 

dy' dzn di/n d^'da da 



sup 



-1 N 

1 + {2'^'(x^ - yn)f + (2fe'(z„ - y„))'J nr=i' [l + 2-^MU,k'r)(^p^^^)M^ 

^ ^2{p+p' + l)]Mi-i)k+il'-^)k' 2l3j{3X ,f^)'kfj^2('^-^)'-'^^^''^^ 

^ (-.2(p+p'+i)j+(i-3)fe+(i'-3)*;'2/3.7 0:fe'.M)-feA'+fe-fc' 
< (j2iP+p'+'^):i+ii-k)k+(i'-k)k'2-(k-k')iJ-i)^ 



where we have used the size estimate for 0i, r > and (6.15) at the first step; the 
expressions (6.21 ) and the support sizes of ajk-, bj'k' at the second step to integrate 



in the order specified; (6.24 1 at the third step; and (6.19) at the final step. Since the 



roles of X and y may be interchanged to get the same bound, the desired conclusion 
is established. □ 



Corollary 6.4. For 5 > \, the statement of Cor. 5.4 holds for A £ IP'' {A, Q 



Proof. The proof is identical to that of Cor. 5.4 



□ 



7. Proof of Theorem 11.11 

Now in a position to prove Thm. [TT| we first treat the classes of operators IP''{A, Cj) 
for j = 0, 2. Since there are standard elliptic 'I'DOs that map any L^-Sobolev space 
isomorphically onto any other, and 



it suffices to assume that s = r = in Thm. 1.1 and show that A g /^''(A, Cj) is 
bounded from ^ L^. These assumptions will be used without further reference 
throughout this section as will the decompositions of the operator A introduced in ^ 
and f|6] The proof consists of combining the estimates for ^oo and A — Ao~ Aoo- 
Of these A^ G \I>p+'^ and hence is a bounded map from i/* — H'^^^ , for all s G M 
and r' > r pp = p + l. Fo r /?' ''(A, Cq), the bounds for Aq and A — Aq — Aoo a re g iven 
in Lemma 

Cor. 16.41 It therefore suffices to show that 



and Cor. 5.4 respectively; for /^''(A, C2) they are in Lemma 6.1 and 



r = min max(ro , ri , Too ) , 



where rQ = p - 



S{1 



is as in (5.4), r^o —p + l and 



+ S{l-l 



\P 



ifl<i 



For ^ > 5, one has ri — Too, so max(ro, ri, Too) — max(p+ ^+S{l + ^),p + l) > p + l, 
with equality being attained if and only ii 5 < {I — ^) /{I + Since the right hand 
side is always strictly larger than | in the range of I being considered, choosing S = 
gives r — p + I in this case. Similarly, if / < — i, then r = max(ro, ri,r^) = p + ^ 
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if and only if (5 > ^/(^ — I)- The right hand side being always smaller than 1/3 in 
the given range of Z, choosing any ^ < S < ^ suffices. For — ^ < I < 5, we get 

1 



min max [ p + SI + 



mm 

(5>i 



mm 

5>i 



-^,P+l + Sil + l),p + l 



p + SI + max 



P + SI + 
l + l 



1 + 35 
6 



1-5 1 + 35 



The case I = is handled similarly. The details are left to the reader. 



□ 



7.1. Remark. Finally, we point out how, for the geometry (A, Ci) arising for the 
normal operator in single source seismic imaging in the presence of fold caustics, 
and more generally for FIOs with folding canonical relations, Thm. 1 1.1 [ follows from 
the proof for (A, Cq). From [32], one knows that the li neari zed forward scattering 
operator F belongs to I^{Csing), where, as described in §^ 1.2 Csing C T*XxT*Y is 
a folding canonical relation, and F*F G /^'°(A, Ci). However, any folding canonical 
relation Cfoui may, by application of canonical transformations xl,XR on. the left 
and right, resp., be microlocally put in Melrose-Taylor normal form [55], so that we 
may assume that Cjoid C Co x Ar.R„-2. In composing C'foid o Cjou, the xl,Xl^ 
cancel, and a calculation shows that 



XR^C^ia ° Cfoid)XR C Xfl'(A U Ci)xR C At.r" U 



X A 



T'l 



with Cq = N*{x2 - 2/2 = - 2/1)^}'; on the operator level, for A £ I"^{Cfoid), 
B e r''{Cfoid), one has B*A e /"+™''"(A, Ci) [32l[l0]. Since XR and its inverse 
correspond to unitary FIOs, L^-Sobolev estimates for /^''(A, Ci) follow from those 

for /P''(A,Cq' X A7-.Rn-2). The coefficient j being irrelevant to the analysis, the 
argument above for the (A, Co) geometry can be repeated. 



8. MiCROLOCAL INVERTIBILTY : A COUNTEREXAMPLE 



In this section we describe an example related to the question of microlocal invert- 
ibility. Given an oper ator A of the type considered in this paper, e.g., 
A G /^''(A,Co) as in (2.6), that is also elliptic in the sense that 

\a{x,y;0)\ 



(8.1) 



inf 

x,y,e 



- > c> 0, 



we would like to determine whether it is possible to left-invert A microlo cally in 
some appropriate 1^ (A, Co) class. More precisely, and in view of Thm. |l.l[ we 
pose the following: 

Question: Let r — r(p,l) denote the regularity exponent obtained in Thm [777] 
Given p, I, does there exist p', I' G M satisfying 

rip,l) + rip',l') = 
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and e > depending o nly o n p,p',l,l' such that for every A e /P''(A,Co) that it 



elliptic in the sense of (8.1 ), one can find B ^ 1^ (A, Cq) such that 



(8.2) E = BA — I maps i?comp boundedly to Hi^"^ for all s G 



We will call an operator E satisfying the Sobolev mapping property in ( 8.2 ) smooth- 
ing of order at least e. The following result shows that the answer to this question 
is no in general. 

Prop osition 8.1. Let % be the Hilbert transform along the cubic in defined in 
{1.1), and let if) ^ {M?) . Then the operator Hq = tpHtp lies in I~2'2 (^A,Cq) but 



there does not exist any B that maps H^^^^p boundedly into H^^^ for every s G M 
and for which BHq ~ I is smoothing of any positive order. 

8.1. Facts about H. Since the operator H is translation-invariant, its Fourier 
transform is a multiplier operator, i.e., there exists a function m on such that 

(8.3) H/(6,6)-"i(ei,6)/(ei,6), with ™(6,6) = /e*'«^*+«^*')|. 

Properties of multipliers such as m have been extensively studied in the literature; 
see [IS], [33] and references there. We note below without proof a few well-known 
facts about m that are crucial to the proof of the proposition. 

(i) m G L°°(R2) n C°°(M2 \ {0}). 

(ii) m(p5i,p3^2) = mill, ^2) for all p > and all ^ = (Ci,6) G I^^- In fact, 
m(/i, 1) is an antiderivative of the classical Airy function composed with an 
afhne transformation. Sec [1; Ch. 10.4] for details. 

(iii) Figure [2] contains the graph of the function fi 1—^ rn{fj,, 1). We observe that the 
graph has a zero, i.e., there exists a ^ such that m(£^i,a^i) = 0. In other 
words, there is a unique cubic curve on which m vanishes. 

Properties (p])-(pli|) imply that there exists a constant C > 



(8.4) |m(C)|<Cmin 



8.2. Construction of the function /q. The key element of the proof of Propo- 
sition 8.1 is the construction of a function /o of fixed compact support which lies 
in and not in H'' for any s > 0, but for which Hf G H^^^ for some e > 0. To 
this end, let us fix x G C{^[—1, 1] such that x > on M, x(0) > 0, set 

Cfe fc-5(i + logfc)-\ fc>l, 

and define 

(8.5) M) = E ^^mi - "fe)x(6 - anl), 

k>l 

where {uk} is a fast growing sequence to be determined in the sequel. In particular, 
we will see that choosing — 2'' works. We collect several facts about /o in the 
following sequence of lemmas. 
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m(^, 1) 




Figure 2. Plot of ^ m{fj,, 1) 



Lemma 8.2. The function fo given by (8.5) has compact support in [— 1, 1] x [—1, 1] 
and lies in L^(M^). Aposteriori, this means that the infinite sum in (8.5) converges 
for almost every ^. 



Proof. Since 



fe>i 



the support property of fo follows easily from that of x- We proceed to estimate 
the norm of /g: 

ll/o|l2 = ll/o|l2= E ^''^fe' f X{^i~^k)m2~anl)x{^,-nk')m2~anl,)d^,d^2 



k,k'>l 



Si + S2, where 



k>l •' 

<C^k-^{l + \ogk)-^ <oo, and 
fe>i 

S2 = E / ^^^1 ~ - anfc)x(6 - nk')x{^2 - anl,) d(, 
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< C ^ CfcCfe' 

k^k' 

< C ^ CfcCfc' 

k^k' 

< C ^ CkCk' 

k^k' 



|?l-rifc|<5|nfe-nfc 



'I -'iCi 



-"fc|>5l"fc-"fc'l 



l?l-"fe'l>jl"fc-"fc'l "'kl-nfcl>5|nfc-nfc'l^ 



I5l-nfcl>3|nfc-nfc'l 



k^k' •' 

<C^CfeCfc.(l + |nfe-nfe/|)-2 

k^k' 

- 1 



(1 + |r7,fe - nfe'|)2 



< C^c\ < oo. 



In the computation above, the second step in the estimation of S2 follows from the 
bound 

sup / x{^2 - anl)x{^2 ~ an\,) d(,2 < C < 00, 

k,k' J 

while the fourth step uses the inequality 1^1 —rifc' I > |n/c~'T'fc'|~ICi~'^fe| ^ 5 1'T-fc— | 
if — rifcl < l\nk — Ti/c'l- This completes the proof of Lemma 8.2 □ 



Lemma 8.3. The function /o given by (8.5) does not lie in Hf^^ for any s > 0. 

Proof. Let V' be any smooth function of compact support. Without loss of gener- 
ality, we may assume that "0 > 0, with ■0(0) > 0. We have to show that -0/0 ^ 
for any s > 0. 

Let us now choose small constants cq, eo > such that 

(8.6) inf{0(O : ICI < eo} > Co, inf{x(ei) : 161 < eo} > cq. 

We compute 



di 



{i + \erfo{v)Mv')dvd7fd^, 



V 



where at the last step we have replaced the domain of {£,,ri,r]') integration by the 
subset 

V:^{{tv,v') ■■ iC-^/l <eo,|e-r/| <eo}, 
and then used the first inequality in (8.6 1. Substituting the expression for /o from 
(8.5) in the last integrand and replacing the infinite sum (of non-negative sum- 
mands) in /c, fc' by the smaller subsum over fc = fc' we obtain 



CkCk' 



k,k'>l 



(1 + l^r)'x('7i - "fe)x(»72 - ant) 



V 
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X xiVi - n-k' )x{V2 - ""-fe' ) dr]' d£, 



k>l 



k>l 



(1+1^1)2 dTjdv'd^, 



where 



Pfc := : le- (nfc,an|)| < so that VkxVkxVkQV, 

and the lower bound for the integrand in the last step follows from the second 



inequality in (8.6). On Vk, the variable <^ satisfies the lower bound 1 + |^| > 
1 + |an|| — 1^2 ^ '^^k\ — 1 ^ ^0 + > so that there exists a constant c > 

for which 



2„6s 
k ■ 



k>l 



The sum on the right hand side above diverges to oo for any s > if rife increases 
sufficiently fast, for instance exponentially. This completes the proof of the Lemma 
lO □ 



Lemma 8.4. There exists Sq > such that Hfo G H'^" . 



Proof. It follows from (|8.3| and (|8.5|) that 

mfowh - 



k,k' •' 

X X(6 - nk')x{^2 - anl) d^ 



Ti +T 



2- 



We now proceed to estimate the two sums separately. For the first, we write 

Ti = E / Hoi'mi - nk)'m - anD^a + d^ 



(8.7) 



k>l 



E- 

k>l 



=:Tii + Ti2, where 



£k ^ ■■ - nk\ < nl \^2-anl\<nl} 
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Here k is a fixed small positive constant (in fact, any k < | will work). The 
multiplier estimate (|8.4[) yields 



Til < 



(8.8) 



But on £k 



k>i 



<C^cl sup 
fe>i ^^^'^ 



6 

<■ 5" 

+ 



2s 



|6 - \ = |6 - a((Ci - rik) + nkf\ 
^ 1^2 - anl\ + \a\\Ci - Ukf 



(8.9) 



while 
(8.10) 
and 

(8.11) 



- 3|a|nfc|^i - n/cP + 3|a|n^|fi - n^l 



<n-, + \a\nl- + 3\a\nl+''^ + 3\a\nl+- 



1 + ICI< 1 + 161 + 161 

< 1 + nfc + lan^l + j^i - n^j + |6 - anW 

< 1 + nfe + I + < Gn\. 



Combining estimates (8.9), (8.101 and (8.111, and using the fact that — n^!^, 
we arrive at the following bound for Tn: 



Til < 



.2„6s V'fc ) „2k 



fc>l 



'fe k>l 



Ci.n 



2 4k+6s-2 



< oo 



if s > is chosen small enough so that 4k + 6s < 2. 
We will now estimate T12. For any iV > 1, 

Ti2 = E^' / H0\'m~nkfm2-anlf{l + \eyd^ 
k>i •'^i 

(i + l6lP(i + l6lPrf^ 



^ 2 /• (1 + ^fc + 16 - »fc|)"'^(l + + 16 - anl\)^^ 
<C^Y.<^l{l^n,r{l + \anl\r [ 



[(l + l6-'^fc|)(l + l6-a?^?.| 

< Cat^ c^n^^-'^^ ({1 + 16 - n,.|)-'(l + 16 - anWy^dS, < 00. 
fc>i •' 



i2Af-2s 



The last step follows from the fact that on £^ at least one of the quantities |6 — "■fel 
or 16 ~ ck'^fel must be > nj;. For the sum in the last step to converge for a given 
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choice of s and k, we must choose N large enough so that 2s — N < —2 and 
8s- kN < 0. 

Finally, we turn to the estimation of T2. This will be done almost exactly in the 
same way as Ti but this time keeping in mind the almost orthogonality of the 
summands. More precisely, 

T2 - ^ ccfe, / |m(oi'(i + lerm - - 

f 



= ^ CkCk' 
k^k' 

k^k' 



ICl-"fel<5|nfc-nfc'l 



ki-"fe'l>3l"fc-"fc'l "'I5l-nfcl>3|nfc-nfc'l. 

2/. . i*i2ns 



k^k' "'l5i-nfc'l>5l"fc-"fc'l 

,3 



< 



CkCk' 



k^k 



; {l + \nk-nk'\y 



X xiCi - nk')x{^2 - ank,)d^ 



We observe that the integrand in the last step above is of the same form as the one 



occurring in the expression for Ti (see the first line in (8.7)); the only distinction is 
that the function in Ti has been replaced by x, which makes no difference to the 
estimation process. We leave the reader to verify that the same set of arguments 
that was used to estimate Ti now yields a number r = min(2 — 4k— 6s, kN—8s) > 
such that 



so that 



To < 



CkCk'Uk 



k^k 



(1 + \nk - Uk' 



< C 



£1 

k^k' 



2 -2r 
Ck^k 



1 + \nk ~ rik' 



£1 

k^k' 



CI, 



1 + |nfe - nk'\Y 



< 00. 



This completes the proof of Lemma [8^ 



□ 



8.3. Proof of Proposition 8.1, It is easy to check that Hq £ /^2'2(A,Co), 
and we leave this to the reader. Without loss of generality, we may assume that 
the function i/j in the statement of the proposition is identically 1 on [—1, 1]^ and 
vanishes outside [—2, 2]^. Let us suppose that there exists an operator B that maps 
^comp boundedly to H(^^ for all s € K and such that BHo = I + E, where _B is a 
smoothing operato r of some order > sq > 0. Then /o = BHo fo — Efo, where /o 
is the function in (8.5 1. Since supp(/o) C [—1,1]^ by Lemma 8.2 it follows from 



Lemma 8.4 that 'H(V'/o) = H/o e H''^ for some small Si > 0, so BUofo & H^^c- 
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On the other hand, Efa G 
fo = BHofo - Efo e 
Lemma 18.31 



So 

loc 

^loc where S2 



since fo E by Lemma 
min(so,Si) > 



This imphes that 
thereby contradicting 
□ 
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